Geodesic scattering by surface deformations of a topological insulator 
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We consider the classical ballistic dynamics of massless electrons on the conducting surface of a 
three-dimensional topological insulator, influenced by random variations of the surface height. By 
solving the geodesic equation and the Boltzmann equation in the limit of shallow deformations, 
we obtain the scattering cross section and the conductivity a, for arbitrary anisotropic dispersion 
relation. At large surface electron densities n this geodesic scattering mechanism (with a oc ^/n) is 
more effective at limiting the surface conductivity than electrostatic potential scattering. 

PACS numbers: 73.23.Ad, 73.25.+i, 73.50.Bk 



I. INTRODUCTION 

Topological insulators such as Bi2Se3 form a new class 
of materials, characterized by an insulating bulk and a 
conducting surfaccjii^ The surface states are massless 
Dirac fermions with spin tied to momentum by spin-orbit 
coupling. Time-reversal symmetry prohibits backscat- 
tering and prevents disorder from localizing the surface 
states. The surface conductivity can therefore be unusu- 
ally large, offering potential applications for electronics. 
The limitations on the conductivity of Dirac fermions 
imposed by random potential fluctuations are well un- 
derstood (mostly from extensive studies of graphene^.). 
Here we study an altogether different non-electrostatic 
scattering mechanism, originating from random surface 
deformations. 

The epitaxial growth of Bi2Se3 films is known to pro- 
duce random variations in the height profile z = C(a;, y) 
of the surfacCi^ These surface deformations correspond 
to terraces of additional layers of the material (of typ- 
ical height H = 2nm and width W = 10 nm). Since 
the Dirac fermions are bound to the surface, they are 
forced to follow its geometry. Like photons in curved 
space-time, the electrons follow the geodesic or shortest 
path between two points, although here the curvature is 
purely spatial)^ (The metric tensor of the surface does 
not couple space to time.) The geodesic motion around 
deformations constitutes a scattering mechanism that by 
its very nature is energy independent, and which there- 
fore is qualitatively different from potential scattering. 

Our problem has no direct analogue in the context 
of graphene. Ripples of a graphene sheet do scat- 
ter the electrons, but this is not geodesic scattering: 
Ripples in graphene are described by gauge fields and 
scalar potentials in a flat space.? Space curvature effects 
may appear around conical defects (pentagon and hep- 
tagon rings), but these are rare in graphenCi^ An early 
study of geodesic scattering in condensed matter that we 
have found in the literature is by Dugaev and Petrov^ 
with possible applications to intercalated layered crys- 
tals. The present work goes beyond their analysis by in- 
cluding the effects of an anisotropic dispersion relation, 
which is a major complication but relevant for topological 
insulators. 



The paper is organized as follows. In Sec.|TT]we inves- 
tigate the classical motion of the surface electrons in the 
presence of surface deformations. The geodesic equation 
is solved in the regime H/W ^ 1 of shallow deforma- 
tions, to obtain the differential scattering cross section 
S. In Sec. mil we use the linearized Boltzmann integral 
equation to compute the conductivity tensor a from S. 
This is a notoriously difficult problem for an anisotropic 
dispersion relation^ In the regime H/W ^ 1 we are able 
to find a closed-form solution, by converting the integral 
equation into a differential equation. Results are given in 
Sec. lIVI and in Sec.|V]we discuss the experimental signa- 
tures that distinguish geodesic scattering from potential 
scattering. 



II. GEODESIC SCATTERING 

A. Geodesic motion 

We consider the surface of a topological insulator in 
the x — y plane, deformed by a locally varying height z — 
C(a;, y). The dispersion relation of a locally flat surface is 
an elliptical hyperboloid. 



E = 



■^xPI + vIpI + vIpI 



(2.1) 



where we have taken the x, y, z axes as the principal axes 
of the elliptical cone. In general, all three velocity compo- 
nents Vx,Vy, Vz may be different. For an isotropic disper- 
sion relation in the x — y plane we have in-plane velocities 
Vx — Vy ~ vf, but the out-of-plane velocity Vz may still 
differ. 

We have included a mass term e in Eq. (|2.ip in order 
to have a nonzero Lagrangian, 



L = ^ XiPi ~ E = -e 1 - y^^{xi/t 



o^ 



(2.2) 



with Xi = dxi/dt — dE/dpi and i = x,y, z. In the final 
equation of motion e will drop out. 

The constraint that the motion follows the surface im- 
plies i = {dC/dx)x + {d(/dy)y, which can be used to 



eliminate i from the Lagrangian. The result can be writ- 
ten in the form 



L 



1 



Vx^g^^i'^i'', 



(2.3) 



with g^y the metric tensor (made dimensionless by 
pulling out a factor w^). Summation over repeated in- 
dices fj-jV — 1,2 = x,y is implied and upper or lower 
indices distinguish contravariant or covariant vectors. 
Explicitly, we find 



2„,2 



ffx. = 1 + (dC/dxYu 
vl, + {dC/dyfv, 



9yy 

9xy 



2„,2 



9yx = idC/dx){dC/dy)vl„ 



(2.4a) 
(2.4b) 

(2.4c) 



where we have abbreviated Vij = Vi/vj. The inverse of 
the tensor g^i,, denoted by g^" , has elements 



gyy^D-'[vl, + {dc/dxfv. 



i/2 

2„,2 

yzi 



g^y = gy- - ~D-'{dc/dx)idc/dyX, 



D^l + (dC/dxfvl, + (dC/dyfu 



(2.5a) 
(2.5b) 
(2.5c) 
(2.5d) 



The Euler-Lagrange equation dL/dx'^ 
{d/dt)dL/dx^^ gives the inhomogeneous geodesic 
equation 



9.10 



■^ ~r -'^ ^11''^ ■^ 



, 1 dL 
T'dt' 



The coefficients F^^^ are the Christoffel symbols. 






\S 



d 



dxt^ 



9su 



dx^^f"- 



(2.6) 



(2.7) 



The nonzero right-hand-side in Eq. (J2.6I) may be elim- 
inated by a reparameterization of time, from i to t such 
that dr/dt = -~L{t)/e. We thus arrive at the homoge- 
neous geodesic equation 



2^A 



d x 



dT^ 



> dx^' dx" 

r^ 

'"^ dr dr 



0. 



(2., 



Since e does not appear in this equation of motion, it 
holds also in the limit of massless electrons. 



B. Scattering angle 

We consider the scattering from a surface deformation 
C(x, y) of characteristic width W and height H large com- 
pared to the Fermi wave length Xp- The scattering may 
then be described by the classical equation of motion, 
which is the geodesic equation (|2.8p . 

An electron with wave vector fc incident on the de- 
formation with impact parameter b at an angle 9^ with 
the X-axis is scattered by an angle 9{0k,b), resulting in 
a differential scattering cross section S{9k,d) = \db/d9\. 




FIG. 1: Geodesic trajectory of an electron deflected by a cir- 
cularly symmetric deformation (characteristic width W). The 
impact parameter b, incident angle Ok, and scattering angle 
9 are indicated. The greyscale background shows the height 
profile of the Gaussian deformation 1)2. lOp . 



Multiple trajectories may lead to the same scattering an- 
gle so that 9{dk, b) cannot be inverted. Then the function 
has to be split into several invertable branches i and the 
cross section becomes S{9k,9) = J2i \dbi{9k,9)/d9\. 

These quantities may be calculated by numerically 
solving the geodesic equation. Analytical progress is pos- 
sible in the physically relevant regime H/W ^ 1 of shal- 
low deformations. As shown in App. IA21 the scattering 
angle is then given by 



^,b)=- 



Ty^{S:,b)di. 



(2.9) 



Here r^^{x,y) is obtained from F^j^(a;, y) by a rotation 
of the coordinate axes over an angle 9k (so that the 
electron is incident parallel to the i-axis). To leading 
order in H/W and b/W the scattering angle scales as 
9 = 0{H'^b/W'^). 

One simple example is the case of a Gaussian deforma- 
tion, 



C(a;,?;) = i/exp[-(a;^ 
which yields (see App. IA3|) 



Z/')/2M^' 



9{9kM 



y±be~^ '^ 



2 W^ 
X {cos^9k+vl^sm^9k), 



(2.10) 



(2.11) 



in the shallow deformation limit. The geometry is de- 
picted in Fig. [T] We will use this example throughout 
the paper to illustrate our general results. 



III. CALCULATION OF THE CONDUCTIVITY 

A. Linearized Boltzmann equation 

We investigate how geodesic scattering influences the 
surface conductivity a of the topological insulator. We 



assume a ^ e^//i, so that we may use a semiclassical 
Boltzmann equation approach. In the presence of an ex- 
ternal electric field E, the occupation fk = fo{Ek)+gk of 
the electron states deviates to first order in E according 
to the linearized Boltzmann equation, 



dfo 
dEk 



e-yfc-£;-^Q(fc,fc')(.9fc-.gfc')- (3-1) 



Here, Vk = d£k/fidk is the velocity and Q{k, k') the 
scattering rate from k to k' (equal to Q{k' , k) because of 
detailed balance). The sum over k' runs over all states of 
the (d-dimensional) momentum space. In the continuum 
limit, J2k ~^ ^ I dk /{2ttY, where V is the d-dimensional 
volume (d = 2 in our case). Spin degrees of freedom do 
not contribute to the sum since the helical surface states 
have definite spin direction. Particle conservation leads 
to the normalization condition 



Efffe^o- 



(3.2) 



the use of q{k, k') — q{k' , k) shows that <t is a symmetric 
tensor. 

For a low density Af of scatterers, the scattering rate 
q{k, k') can be related to the differential cross section S 
of a single scatterer (averaged over all scatterers). In the 
two-dimensional case of interest here, the relation is 



N\vk\S{6k,6k')d6k' ^ q{k,k') ^ '^^'P 



{2T:Y\hvk'V 



(3.9) 



where Ok is the angle between Vk and the x-axis. The 
Eq. (j3.6p which determines the vector mean free path 
then takes the form 



J^\vk\ I d0k'S{0k,9k'){Ak- Ak')^Vk. 
Jo 



(3.10) 



For the solution of this equation (and the interpreta- 
tion of the results), it is convenient to follow Ziman^iii^ 
and define an anisotropic relaxation time T(fc) by 



The electric field can be eliminated from Eq. p.ip by 
means of the vector mean free path A^ , defined by^i^ 



gk 



dfo 



eE ■ Afc, 



dEk 
Y,Q{k,k'){Ak-Ak' 



= Vk- 



(3.3) 

(3.4) 



For elastic scattering, Q{k, k') = S{£k^£k')q{k, k'). Us- 
ing dk — dk^dSp — dSkdSp/lfiVkl, with dSp a Fermi 
surface element, Eq. p.4p can be rewritten in terms of 
the density of states N{£jr) at the Fermi energy. 



N{£f) = (2^) 



-d 



dSp Ihvk 



(3.5) 



The integral ^ dSp extends over the Fermi surface. The 
result is 

VN{£F){q{k, k'){Ak - Ak'))k' = Vk, (3.6) 

with (• • • )fc denoting the weighted average over the Fermi 
surface. 



{f{k))k = 



' dSp fik)\hvk\~^ 
SdSF\hvk\-^ ■ 



(3.7) 



The normalization condition (|3.2p becomes {Ak)k = 0. 

At zero temperature, the conductivity tensor is given 
by 



o- = — 2_] S{£k - £f) Vk ® Afe 

k 

= e^N{£F){vk®Ak)k. 



(3.8) 



The direct product ® indicates the dyadic tensor with 
elements [t)fc]i[Afc]j. Substitution of Eq. p.6p for Vk and 



'{k) 



VN{£f){{1 - Vk ■ Vk')q{k,k'))k'. (3.11) 



Using Eq. (13.91) this can be rewritten as 

<.27r 

= KfUuA I 

T{k) 



kT=^\vk\ f d9k'Si0k,ek')[i-cos{ek'-ek)]. 

.'^) Jo 

(3.12) 



B. Isotropic dispersion relation 

For isotropic dispersion relations (when Ek depends 
only on |fc|, so that the velocity v = VFk is aligned with 
the wave vector), the linearized Boltzmann equation can 
be solved exactly^ This applies, for example, to surfaces 
perpendicular to the [111] direction of Bi2Se3. We con- 
sider this simplest case first. 

Since the deformations do not have a preferred orienta- 
tion and the dispersion is isotropic, the average scattering 
cross section S{9k,0k') only depends on the scattering 
angle = 6k — Ok' , independently of the incident direc- 
tion. The solution to Eq. p.6p is then Afc = rvk with a 
relaxation time r given by 

1 Z'^'" 

-:^Uvf d0S(9)(l'~cose). (3.13) 

T Jo 



Substitution into Eq. (|3.8p leads to a scalar conductivity 
a given by the Drude formula, 

...^„„,,4l4|l. (,I4) 

In the second equality we inserted the density of states 
N{£f) = £f/{'^t^^^v\) of a Dirac cone with a circular 
cross section. 



The regime H/W <C 1 of shallow surface deforma- 
tions is characterized by predominantly forward scatter- 
ing {\9\ <C 1). Then the relaxation time (|3.13|) is given 
by the second moment of the scattering angle, 



- ^^J\fvF f des{0)e^ 



(3.15) 



We substitute the relation S{9) = {\d9{b)/db\~^), where 
(• • • ) indicates an average over the (randomly oriented) 
scatterers. The integration over scattering angles 6 be- 
comes an integration over impact parameters b, 



- = WVF 

T ^ 



dbe^ib) 



(3.16) 



From Eq. ()2.9p we infer the scaling 1/t (xWx [H/W)^ 
of the relaxation rate with the characteristic height and 
width of the surface deformations. (The additional factor 
of W comes from the integral over b.) This scaling was 
first obtained by Dugaev and PetrovJ. Eq. (I3.14p then 
gives the scaling of the conductivity 



<T = constant x 



e'Ep 



1 W^ 



h n AfvF H^ 



(3.17) 



C. Anisotropic dispersion relation 

We now turn to the case of an anisotropic dispersion 
relation. There is then, in general, no closed-form solu- 
tion of the linearized Boltzmann equation. ^^ One widely 
used approximation for the conductivity, due to Ziman,— 
has the form 



CZiii 



e^N{£F){vh<^VkT{k))k, 



(3.18) 



with T{k) the anisotropic relaxation time p. lip . As we 
will show in the following, this is a poor approximation 
for our problem, but fortunately it is not needed: In 
the relevant limit H/W ^ 1 of scattering from shallow 
surface deformations an exact solution becomes possible. 
For shallow deformations forward scattering dominates, 
1^1 — \0k — ()k'\ ^ 1- This allows for an expansion of A^/ 
around 6*^, which reduces the integral equation p.6p to 
a differential equation. 
With the notation 



Afp(0) 



des{(t),(t) + e)ep, (3.19) 



the expansion to second order of Eq. p.lOp can be written 
as 

Mi(0)^A(0) + \M^{cp)^\{cp) = -^e^t (3.20) 

We introduced a complex variable A = Aj, + iKy to com- 
bine the two components of the vector mean free path. 
Denoting the radius of curvature of the Fermi surface by 



k((/)) = dSp/dcj), the normalization condition (13.21) be- 
comes 



27T 






A((/.) = 0. 



(3.21) 



Once we have the solution of Eq. p.20p , the conductivity 
tensor elements follow from 



2 /-^TT 7 i 

fT,, ± ayy = ^Re / -^ e^'^K{^)X{cl>), (3.22a) 
o.y ^ay^^^-^lmT ^ e'^K{c^)X{(f>). (3.22b) 

A further simplification is possible if the average scat- 
tering angle vanishes, Mi{ip) = 0. Then the second mo- 
ment M2{4>) of the scattering angle is, within the for- 
ward scattering approximation, directly related to the 
anisotropic relaxation time: 



1 



r(0) 



= Wvicj))M2{^). 



(3.23) 



Eq. (|3.20p can now be solved in terms of the Fourier 
transforms 



2^ 

2-!T 



■e-""^w(0)T(0), (3.24a) 

K„ = r ? e-™'^^(0), (3.24b) 
Jo 27r 

Xn^ j'^ '^e-"^^\{<t>), (3.24c) 



resulting in 



A„ = — - — h constant x (5„ o- 



(3.25) 



The normalization constant can be determined from Eq. 

(EUD. 

Inserting the solution into Eq. p.22p we obtain the 
conductivity 



Gxx ± (yyy = — Re } ^ 



n— — oo 



^ xy — ^yx 



ii- E 



-l^-n±l 
9 ■ 



^_lK_„_l 



(3.26a) 
(3.26b) 



n— — oo 



For simplicity we have assumed an inversion symmetric 
Fermi surface, for which k±i = so that the normal- 
ization constant in Eq. (|3.25p does not contribute to the 
conductivity. 

In the case of an isotropic Fermi surface, only the 
Fourier components ^o = vpT and kq = kp arc nonzero. 
From Eq. (|3.26p . we then find Uxy = = Gyx, (Jxx — 
ayy = {e'^/2h)kpVFT, in agreement with Eq. p. 141) . 



Comparing with the Ziman approximation p. 181) for 
the conductivity in terms of the anisotropic relaxation 



time, we see that it can be written in the same form 
p.26p . but without the factor l/n?. It therefore deviates 
strongly from our forward-scattering hmit, except in the 
case of an isotropic Fermi surface (when only n — 1 con- 
tributes). 



IV. RESULTS 

A. Isotropic dispersion relation 

In the shallow deformation limit the conductivity is 
given by Eq. p. 171) , up to a numerical prefactor of order 
unity. We have calculated this prefactor for Gaussian 
deformations of the form (I2.10p . randomly distributed 
over the surface. We assume that the deformations are 
shallow, H/W <C 1. For simplicity, we also take the 
same parameters H and W for each deformation. From 
Eqs. dmH), (jXTil) . and ((XTCl) we obtain the resuh 



16%/2 Ef 



W^ 



7r0r hvpM [Hvp/vzY h 



(4.1) 



The factor vf/vz is there to allow for an out-of-plane 
velocity Vz that is different from the in-plane velocity 
Vx = Vy — vp- The result (j4.ip confirms the scaling 
behavior (|3.17|) and gives the numerical prefactor. 

To relax the assumption H/W ^ 1 of shallow defor- 
mations, we solved the geodesic equation (12. 8p numeri- 
cally for the Gaussian case. The corresponding Christof- 
fel symbols were taken from Eq. (IA2|) with Vx = Vy — vp- 
Using the scattering angle 9{b) that we obtained from the 
numerics, we calculated the conductivity following from 
Eqs. (l3l3l[3l4D . 

As shown in Fig. [21 the numerical results deviate from 
the scaling (|4.ip only for relatively large ratios H/W > 
0.5. The deviations are oscillatory, due to electron tra- 
jectories that circle around the deformation as depicted 
in the inset (b) of Fig. [51 Inset (a) shows generic tra- 
jectories for electrons scattering off a shallow Gaussian 
deformation. Notice the focussing of trajectories as an 
analogue of gravitational lensing. 



B. Anisotropic dispersion relation 

As an example of an anisotropic dispersion relation, 
we consider elliptic equi-energy contours £k = h{v'^k'^ -^ 



Vyky)^''^ with principal axes x and y. As in the previ- 



ous subsection, we investigate shallow Gaussian surface 
deformations. These have zero average scattering angle, 
Mi{(p) — 0, and second moment 



The coefhcient C is given by 
16^/2 



^Ix COS^ . 



c 



W^ 



Hi H-^vyv^ ' 



(4.2) 



(4.3) 



Us 10" - 




H [nm] 

FIG. 2: Surface conductivity of a topological insulator as a 
function of the height H of randomly positioned Gaussian 
deformations (width W = 10 nm, density N = 0.1 W~'^). We 
took an isotropic dispersion relation, with in-plane velocities 
Vx = Vy = vf = 5-10^ ni/s, and a smaller out-of-plane velocity 
Vz — vf/3. The Fermi energy is fixed at £f = 150 meV. As 
discussed in Sec. [Vl these are realistic parameter values for 
the [111] surface of Bi2Se3. Dots represent numerical results 
whereas the line shows the shallow deformation limit (14.11). 



From Eq. (|A12|) we deduce that Eq. (|i^ actually holds 
more generally for any circularly symmetric deformation, 
the only difference being in the expression for C. 

Using Eqs. (|3.23l) and (I3.24ap one obtains the Fourier 
coefficients 



fe.^g(;^^V""''^+i"iv+.;.) (,„ 



Af \l + v„ 



for n even, and zero for n odd. The elliptic dispersion 
relation leads to 



«(0) 



£f 



hvx (sin (j) + vl^ cos^ 0)3/2 



(4.5) 



The Fourier coefficients k„ are also nonzero only for n 
even. (Since their expressions are rather lengthy, we do 
not list them here.) 

From Eq. p.26p we find that the off-diagonal compo- 
nents of the conductivity tensor vanish, while the diago- 
nal components are given by 

^{"} ^ T ZJ 0~2(^"+l ±^n-l)(K„+l ± Kri-l). (4.6) 
n>l 

The series converges rapidly. 

The ratio crxx/f^yy depends only on the anisotropy 
Vyx = Vy/vx- It is plotted in Fig. [31 For comparison, 
we also show the Ziman approximation fJ^-^^^^^^ (obtained 
from the forward-scattering limit (|4.6p without the l/n^ 
factor). As expected, it deviates substantially upon in- 
creasing the anisotropy (notice the logarithmic scale). 




VyjV^ 



FIG. 3: The solid line shows the ratio of conductivities 
i^xx/o-yy as a function of anisotropy Vy/vx, calculated from 
Eq. (|4.6p . The dashed line corresponds to the Ziman approx- 
imation. 



V. COMPARISON WITH POTENTIAL 
SCATTERING 

A. Carrier density dependence 

The energy independence of the mean free path I — 
vpT is the hallmark of geodesic scattering. It implies the 
square root dependence dependence a ex ^/n of the con- 
ductivity on the surface electron density n. This follovifs 
from Eq. p.l7p vifith Ep = hvp^^T^n for an isotropic 
Dirac cone, or more generally from the scaling a ex Sp 
for a noncircular Fermi surface (of area Sp oc ^Jn) . 

As discussed in the context of graphene,— i^^ electro- 
static potential scattering typically gives a faster in- 
crease of the conductivity with increasing carrier density. 
Coulomb scattering from charged impurities and reso- 
nant scattering from short-range impurities both give a 
linear increase a on n (up to logarithmic factors). Scatter- 
ing from a potential landscape with a Gaussian correla- 
tor gives an even more rapid increase a ex n'^". Geodesic 
scattering, with a ex n^/^, would therefore form the dom- 
inant conduction-limiting scattering mechanism at high 
carrier densities. 

For a quantitative comparison of geodesic and poten- 
tial scattering, we consider the [111] surface of Bi2Se3 
with Gaussian deformations given by Eq. (|4.ip . We take 
isotropic in-plane velocities Vx = Vy = vp = h ■ 10^ m/s 
and a smaller out-of-plane velocity Vz = vp/i}^'^^ We 
adopt the following numerical parameters for the de- 
formations from an experimental image— characteristic 
width W — 1{) nm and height H — 2 nm, covering 40% 
of the surface area so A/" = 10^^ cm~^. The carrier den- 
sity dependence of the conductivity for geodesic scatter- 
ing, following from Eq. (|4.ip . is plotted in Fig. 0] (solid 
curve) . 

To compare the geodesic scattering to typical poten- 
tial scatterers, we also show the corresponding results for 
scattering from charged impurities (dashed) and Gaus- 
sian potential fluctuations (dotted) in Fig. 2] 
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FIG. 4: Conductivity as a function of carrier density. The 
influence of three different sources of scattering is shown: sur- 
face deformations (solid line), unscreened Coulomb impurities 
(dashed line) and Gaussian correlated potential fluctuations 
(dotted line). The parameters used for the plot are given in 
the text. 



For charged impurities (charge Q = e) we con- 
sidered the unscreened Coulomb potential U{r) = 
{Qe/4:TTef)er)\r\~^, as the extreme case of a long-ranged 
potential. We took e^ = 80 as a typical value for the 
dielectric constant and kept the other parameter values 
as before. The semiclassical conductivity is then given 
by^^ 



e^ n 27rh?Vp 



h K 



2 , "0 
c "0 



Qe 



(5.1) 



For Fig. m we used Afc — 2.5 x 10^-^ cm^^ as the density 
of impurities. 

For a potential landscape with Gaussian correlator 
(range ^, dimensionless strength Uq), 



{U{r)U{r')) 



U^{hvpf 



■ exp 



r — r 



2<2 ^ V 2^2 

the conductivity takes the functional fornJ^ 

_ e2 47r<2g47r«C^ 



(5.2) 



(5.3) 



(The function /i is a Bessel function.) For Fig.|4]we took 
Uo = 0.1 a.nd£, = W = 10 nm. 

The parameter values used in Fig. |4] are only for the 
purpose of illustration, but the point to make is that 
geodesic scattering dominates over potential scattering 
for large carrier densities. 



B. Anisotropy dependence of conductivity 

In the case of an anisotropic (elliptical) dispersion re- 
lation the conductivity will be direction dependent. This 
situation arises for example if the surface of Bi2Se3 is not 
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FIG. 5: Ratios of conductivities along the two main axes of 
the dispersion relation are shown as a function of anisotropy 
Vy/vx- The influence of four different sources of scattering is 
shown: surface deformations (solid line), unscreened Coulomb 
impurities (dashed line), Gaussian potential fluctuations (dot- 
ted line), and short-ranged potentials (dot-dashed line). The 
parameters used for the plot are given in the text. 



^ — >■ (at fixed n) of a short-ranged potential. 

From the double-logarithmic plot in Fig. [S] one can 
see that there is an approximate power law dependence, 
(^xx/f^yy oc {vy/vx)~^, ovcr at least one decade. The 
exponent is p « 3.3 for geodesic scattering, while p = 
2 for short-range potential scattering. Scattering from 
long-ranged Coulomb impurities or from medium-ranged 
Gaussian potential fluctuations gives p <2. 

Anisotropic charge transport in the presence of un- 
screened Coulomb impurities for an elliptic dispersion 
relation was also discussed in the context of strained 
graphenei^ There it was argued that (Jxx/o'yy oc 
{vy/vx)~'^ on the basis of a power-counting argument. 
Our numerical solution of the Boltzmann equation gives 
a smaller exponent p sa 1.3 in that case. 

To conclude, charge transport dominated by surface 
deformations has a much stronger anisotropy dependence 
than that governed by impurity potentials. This highly 
anisotropic transport behavior is a distinct characteristic 
of geodesic scattering. 



in the [111] direction. Geodesic scattering implies a cer- 
tain universality for the directionality dependence of the 
conductivity, if we may assume that the surface deforma- 
tions are shallow {H/W <^ 1) and without a preferential 
orientation (circularly symmetric on average). The ratio 
^xxl^^yy is then only a function of Vy/vx, independent 
of other parameters (such as electron density or density 
and height of the deformations). This universal function 
is plotted in Fig. [3] (solid curve). 

In Fig. [5] we compare this result for geodesic scat- 
tering with corresponding results for potential scatter- 
ing. Three typical impurity potentials are considered, 
of different range: long-ranged unscreened Coulomb po- 
tentials, medium-ranged Gaussian potential fluctuations, 
and short-ranged potentials. The conductivities are ob- 
tained following the general approach of Ref. [I^, by first 
computing the transition rates in Born approximation 
and then solving numerically the linearized Boltzmann 
equation. We took the same material parameter values 
as in the previous subsection. 

The unscreened Coulomb potential gives a ratio 
(Jxx/'^yy which depends only on Vy/vx (dashed line). For 
Gaussian potential fluctuations, the ratio Oxx/c^yy is a 
function of both Vy/vx and n. It is plotted as a dotted 
line in Fig. [5] for n^^ = 1. (If ^ = VF = 10 nm this 



corresponds to the carrier density n 



10 



12 
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Appendix A: Calculation of the scattering cross 
section 

1. Christoffel symbols in rotated basis 

In order to calculate the scattering angle in the geome- 
try of Fig. [2 it is convenient to rotate the coordinate axis 
in the x — y plane such that the electron is incident par- 
allel to the X-axis. Under the linear transformation from 
x,y to X = X cos 0k + y sin 9k, y = —x sin 9^ + y cos 9^ , 
the Christoffel symbol F^^ transforms to 



f^.(i,y) 



dx y dx^' dx" 



dx^' 



dx>^ dx'' 



(Al) 



the same figure we also plot (dot-dashed line) the limit 



Using the expressions (12.41) , (I2.5P , (|2.7p for metric ten- 
sor and Christoffel symbols, we arrive at 



J^X 
fJ,U 



D- 



d\ 



dxf^dx'' 






^lz^-{-"l,~vl^)sui9, 



dC 



V r)-r, 



dy 



(^'L-^'L)sin6'fc 



dC . 



fc I ^ sm t/fc 



77:7 cos 9k 
ox 



■^ cos 9k 
ay 

-^ sm ( 
dy 



(A2a) 
(A2b) 



The factor D from Eq. ()2.5dl) . written in terms of the rotated coordinates, reads 

2 



D = l 



-;— COS 9k- TTZ sm ( 

ox oy 



9C . n dC ^ 
-^ sm 6'fc + — r cos 0k 
ox oy 



(A3) 



r 



The Christoffel symbols (IA2[) appear in the geodesic 
equation for the rotated coordinates, 



2~X 



d X 



dT^ 



r 



dxi" dx'' 
^^ dr dr 



0. 



(A4) 



2. Geodesic equation for shallow deformation 



The geodesic equation (IA4p can be considerably sim- 
plified in the shallow deformation limit H/W ^ 1. Let 
us consider a particle incident on a deformation along 
the i-direction from — oo with impact parameter b and 
velocity 



VxVy{v'lcos^ 



9k +vlsm'^ 9k) 



'1/2 



(A5) 



Since the derivative dy/dr is smaller than dx/dr by a 
factor {H/WY, we can drop this derivative from the 
geodesic equation. The result is 





d^i 
dr^ 


^ '^ XX 


[d^J 


= 0, 


(A6a) 




<Py 
dr^ 


' '- XX 


(diV 
ydr) 


= 0. 


(A6b) 


Furthermore, 


since i 


ix/dr 


= v[l + 


0{H/W)% 


we can 


write d/dr — 


vd/dx. 


This leads to 










dx"^ 


XX' 




(A7) 


The scattering an 


gle 9 


< 1 is 


obtained from 9 — 


lim2_^oo dy/d: 


i, hence 









9fc,6) = - 



^ XX "'•^ 



Inserting Eq. (jA2bl) into Eq. (|M)) and noting that D = 
1 + 0{H /WY , we obtain the scattering angle to leading 
order in H/W, 



e{9k,h) 



We abbreviated 



dx 



dy 



dC\ d^C 



dx I dx^ 



(A9) 



y^b 



a — Vy^ cos^ 6k + v1^ sir? Ok, 
7= (wL-^'L)sin6lfcC0s6»fc. 



(AlOa) 
(AlOb) 



3. Circularly symmetric deformation 

For a circularly symmetric height profile C(x,?/), de- 
pendent only on r = -y/x^ + j/^ — ^J ff- -|- y"^ ^ the term 



proportional to 7 in Eq. (|A9[) vanishes (because it is an 
integral over an odd function of x). The expression for 
the scattering angle thus simplifies further to 



9{9k,b) 



dx 



dy dx"^ 



(All) 



For the Gaussian deformation (I2.10p we obtain the scat- 
tering angle (|2.1ip given in the main text. 

The entire dependence of the scattering angle 9 on the 
angle of incidence 9k is contained in the prefactor a. This 
implies that the moments Mp = j dh9'P of the scattering 
angle depend on the angle of incidence as 



Mp{9k) 



Cr,(T 



Cpi;P,(sin2 0fc-ft;2^cos2 0fc)P, (A12) 



( A8) with Cp a coefficient independent of 9k ■ 
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